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We analyze the vortex dynamics and vortex-vortex interactions in Bose-Einstein condensates con- 
fined in toroidal traps. We show that this particular geometry strongly distorts the vortex dynamics. 
The numerically calculated vortex trajectories are well explained by an analytical calculation based 
on image method and conformal mapping. Finally, the dissipation effects are discussed. 
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Introduction 

The realization of Bose-Einstein Condensation (BEC) 
in trapped dilute alkali gases [jl], ||, ||, [|, ||, ^| has opened 
the possibility to explore fundamental effects of low tem- 
perature physics in a very controllable way. In this sense, 
several effects related to superfluidity have been re- 
cently reported || |) . Among the phenomena related to 
superfluidity, the possibility of observing vortices with 
quantized circulation |lC| ] has aroused a big interest. 
Due to the diluteness of the atomic condensates, accu- 
rate theoretical predictions are possible. In addition, in 
these systems, the diameter of the vortex core is typ- 
ically three orders of magnitude larger than for super- 
fluid Helium, allowing for a more straightforward optical 
observation. Numerous theoretical schemes have been 
suggested for the generation of vortices in BEC, see e.g. 
Refs.Jn], §@g|,0,0. So far only two tech- 
niques have successfully generated vortices in conden- 
sates, namely the stirring of a condensate with a blue 
detuned laser jl^, [T^] , and the phase imprinting method 
achieved via the coherent interconversion between the 
two components of a binary condensate p0[ . In addi- 
tion the creation of a vortex ring after the decay of a 
dark soliton has been recently reported pl| . The dy- 
namics and stability of vortices in trapped Bose-Einstein 
condensates has recently also been a subject of active in- 
vestigation n n, n . 

In the majority of current BEC experiments, the con- 
densate is created in a cylindrical trap, either pancake- 
shaped, or cigar-shaped. However, in the recent years 
the possibility of creating a BEC in toroidal traps has 
aroused a growing interest. Toroidal traps have already 
been used in early Sodium BEC experiments at MIT [g| , 
and in recent experiments on vortices at JILA j2(J . Also, 
very recently the idea of creating a toroidal optical dipole 
trap using Laguerre-Gaussian beams has been proposed 
p6[ . The physics of condensates in toroidal traps have 
been extensively investigated during the last recent years, 
including the ground state properties and elementary ex- 
citations [^t], £8| , the possibility to obtain vortices which 
due to the toroidal topology could be metastable beyond 
some critical nonlinearity |27], |28|, ^| , the generation 
and properties of solitons and the so called "svortices" 
in this kind of traps ||3l| , and even the possibility to 
obtain a mode-locked atom laser with toroidal geome- 
try fl32|. Very recently, we have proposed that a time- 



averaged potential produced by the shaking of an usual 
Gaussian trap also leads to the formation of a double 
well or a toroidal condensate Q. Additionally, if dur- 
ing the condensate shaking a proper phase is imprinted, 
dark solitons and eventually vortices can be created in 
a controllable way within the torus. In Ref. |33|, it was 
already reported, that the particular toroidal geometry 
significantly modifies the vortex dynamics, as well as the 
vortex-vortex interaction. It is the aim of this paper to 
extend the analysis of Ref. Q. We analyze in detail 
the effects of the geometry for toroidal box potentials, as 
well as those produced by the smooth transversal confine- 
ment. Interestingly, we show that the vortex dynamics 
obtained from direct numerical simulations can be ana- 
lytically well described by means of the so-called image 
method, and conformal mapping techniques. Finally we 
analyze the vortex dynamics in presence of dissipation. 

The structure of the paper is as follows. In Sec. I we 
describe the physical system, as well as the numerical 
techniques employed. In Sec. || we consider the case of 
a toroidal box potential, and introduce the image tech- 
nique, as well as the corresponding conformal mapping, 
to analytically describe the vortex dynamics, as well as 
the vortex interaction. The case of a toroidal trap with 
harmonic transversal confinement is considered in Sec. 



Ill, where the effects caused by the smooth trapping 



potential are considered. Sec. 



IV is devoted to a brief 
discussion of the dissipation effects in the system. We 
finalize with some conclusions in Sec. |v|. 



I. MODEL 

At sufficiently low temperatures, the dynamics of 
a condensate in an external potential V can be well 
described by the so-called Gross-Pitaevskii equation 



ih— = -—v 2 * 

dt 2m 

4irh 2 aN 



(1) 



where g — 47r ^ Tt aJV denotes the coupling constant, m 
refers to the atomic mass, a is the corresponding s- 
wave scattering length, and N is the number of con- 
densed atoms. We shall consider in the present paper 
trapping potentials in which the condensate is strongly 
confined along one direction, so that the mean field en- 
ergy is smaller than the excitation energy of the trap in 
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this direction. In this situation, the BEC can be con- 
sidered as quasi-2D, and we can restrict ourselves to the 
analysis of the corresponding 2D GPE, with a modified 
coupling constant g^D — 9 J dz\ip(z)\ 4 , where ip{ z ) is the 
wavefunction in the strongly confined z-direction. Both 
2D |3 and quasi-2D |§ |§ BEC has been recently 
observed in experiments. The potential V in the remain- 
ing two dimensions has an annular shape, which we shall 
consider possessing rigid box boundaries in Sec. ||, or 
harmonic transversal confinement in Sec. III. We restrict 



ourselves in the following to the situations in which the 
BEC can be considered in the Thomas-Fermi regime, 
i.e. the mean-field energy is much larger than the typ- 
ical energy spacing of the annular trap. We consider in 
the following numerical simulations the case oi N = I 0° 
condensed Sodium atoms (a = 2.75 nm). 

We perform the integration of equation (|f|) by an alter- 
nating direction implicit (ADI) Crank-Nicholson scheme 
po! , using cylindrical coordinates, due to the particu- 
lar symmetry of the trapping potential. In the present 
paper, we shall not consider the issue of the experimen- 
tal creation of the vortices, but rather we shall concen- 
trate on the vortex dynamics once created. In order to 
numerically create the vortices, we employ the Crank- 
Nicholson procedure in imaginary time, while ensuring 
at each time step the proper phase pattern necessary to 
create the vortices j4j|. This leads to the creation of 
vortex states, which are free of any other kind of exci- 
tation, and, therefore, allow for a clear analysis of the 
vortex dynamics. The proper phase pattern imprinted 
during the imaginary time evolution has to be carefully 
chosen. In particular, the boundary effects have to be 
taken into account in the method, by adding the corre- 
sponding phase contribution of the image vortices (see 
next section). Once the vortices are created we perform 
a real time evolution. 



II. BOX-LIKE TOROIDAL TRAP 



In this section, we investigate the evolution of vortex 
states in an annular trap with rigid boundaries. In such 
a trap the condensate density in the bulk is constant and 
shrinks to zero at the condensate boundaries on a length 
scale of the order of the healing length 77 = 1/ ^/8nna , 
where n is the condensate density. Due to the absence of 
the potential in the interior of the torus and the homo- 
geneity of the density inside the box, the vortex dynam- 
ics is completely determined by the toroidal geometry. 
In our numerical simulations we have considered a trap 
with inner rigid boundary at a radius R\ = 12/zm, and an 
outer one at R2 — 18/im. We shall not consider, neither 
in this section nor in the next one, dissipation effects. 



Single vortex dynamics 

We consider first the situation in which a single vortex 
is placed in a toroidal trapped condensate. Due to the 
boundary effects the vortex will not remain at rest, but 
on the contrary it will move at a fixed radius around the 
torus, where the direction and the modulus of the vortex 
velocity depends on its circulation k, and on the radial 
position at which it is created. In Fig. El we present the 
numerical results for the vortex velocity as a function of 
the radial vortex position. 

Such dynamics can be well understood by employing 
the so-called image method |^3|, ^4], ^5). Let us briefly 
review this method. In a homogeneous(untrapped) con- 
densate the phase of a state containing a q times charged 
vortex is given by e lqLp , where ip is the polar angle cen- 
tered at the position of the vortex [ff2j . Therefore, the 
superfluid velocity becomes vsf = hVS/m — Tiq/mp e v , 
where S is the condensate phase and p is the distance 
from the vortex line. In the following, we shall only dis- 
cuss vortices with one quantum of circulation since, in the 
cases considered, vortices with q > 1 will decay into q = 1 
vortices. For trapped condensates in the Thomas-Fermi 
regime, the hydrodynamic form of the Gross-Pitaevskii- 
equation allows to approximately impose for the super- 
fluid velocity field the boundary condition of ordinary 
hydrodynamics, namely the vanishing of the normal com- 
ponent of the superfluid velocity field at the condensate 
boundary E| . This boundary condition and the speci- 
fication of the positions and circulations of the vortices, 
leads to the determination of the superfluid velocity field 
p3. Moreover, the vortex dynamics in the homogeneous, 
and also in the boxlike trapped quantum gas, reproduces 
the classical vortex dynamics, which is governed by the 
Kelvin theorem, i.e. a vortex moves with the same veloc- 
ity than the superfluid at its position |Q . This velocity 
can be calculated by employing a well-known method 
in classical hydrodynamics, namely the replacement of 
the boundaries by additional fictitious vortices, so called 
image vortices. These vortices must be placed at the 
appropriate positions with the appropriate circulations, 
such that the resulting velocity field obeys the boundary 
conditions for the considered geometry. 

For the toroidal geometry in question, we have to take 
into account two (in principle infinite) families of im- 
age vortices. For the first family, an image vortex is 

placed at ri m i — -^r^ with circulation —k and one im- 
age vortex with circulation k at the origin |iq| , where 
r denotes the position vector of the real vortex in a 
frame centered at the origin of the torus 46 . The im- 



age vortex at Ti m \ induces a velocity component nor- 
mal on the outer circle, so in addition we need an im- 



These effects will be discussed in Sec. IV 



„+r with circulation n. 

.1 r Rf 

In turn this image will induce some normal component 
on the inner circle, so two more image vortices are re- 
quired: one at r im3 = — —- = with circulation 

— k the other one at the origin with circulation k and 



age vortex at 
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so on and so forth. For the second family, one consid- 
ers an image vortex at r'i m i = with circulation 
—K. This image vortex gives a normal component of the 
velocity on the inner circle so we need in addition one 



r with circulation 



image vortex at r'i m 2 

k as well as one image vortex at the origin with circu- 
lation —k. For this we require another image vortex at 



T im3 



-T?r-- with circulation 

R^r r 



-k and so on. 



The total superfluid velocity is made up by the contri- 
butions of all vortices (both real and image): 



VSF 



\ - n _ r — r n 
,Qn — e z x — — —j. 

m \r — r n \ z 



(2) 



In absence of friction, the superfluid velocity at the vortex 
position equals the velocity of the vortex line. Therefore 
we get for the vortex velocity: 



Vortex-vortex interaction 



We can at this point analyze the vortex-vortex inter- 
action within an annular trap with rigid boundaries. In 
a homogeneous condensate the vortex-vortex interaction 
depends on their relative circulation. In absence of dis- 
sipation, vortices with the same circulation orbit at a 
fixed distance from the central point between them. On 
the other hand, vortices with opposite circulation move 
parallel to each other forming a so-called vortex pair, 
which moves as a whole with a velocity which depends 
on the separation between the vortices. Two vortices in a 
toroidal box exhibit a completely different behavior. We 
shall consider in the following several possible scenarios. 
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where R\ «-> R2 denotes the same term as the precedent 
one, but with Ri and R2 interchanged. This velocity 
field is in good agreement with that calculated from our 
numerical simulations. 
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Vortices with the same circulation created at different 
radii but at the same polar angle show a dynamics very 
similar to that of the homogeneous case: they perform de- 
formed orbits around each other but with an overall drift 
around the torus (see. Fig. ||). On the contrary, vortices 
with opposite circulation created on the same radial co- 
ordinate, but at opposite positions in the torus show an 
interesting evolution, as already discussed in Ref. Q. 
As shown in Fig. ||, the vortices move towards each other 
keeping the radial coordinate constant. When approach- 
ing, they move radially, and separate on a new constant 
radius, i.e. effectively they repel each other. When they 
meet at the other side of the torus, they repel each other 
again, and return to the original radial coordinate, i.e. 
the overall trajectory for each vortex is closed. If the 
vortex-anti vortex-pair is created such that the vortices 
are close to each other at the same polar angle, one vor- 
tex moves very quickly (during a few ms) towards the 
condensate boundary, where it is destroyed. The other 
one moves slightly, changes its initial radial position, and 
later on moves on a fixed radius (Fig. ||). This effect ex- 
plains the results of Ref. p3[ , where a dark-soliton cre- 
ated in a toroidal geometry seemed to decay into just one 
vortex, apparently violating the conservation of the an- 
gular momentum. The dark-soliton actually decays into 
two vortices, but one of them is rapidly destroyed at the 
condensate boundary, while the remaining one is set into 
an orbit around the annulus. 



FIG. 1: Velocity of the vortex in a box like toroidal trap 
as a function of the radial coordinate. The three curves are 
obtained by direct evaluation of Eq. (^) (dashed), Eq. ( |2l| ) 
(dashed-dotted) and by solving the GPE numerically (solid). 



Employing, as above, the image method technique, the 
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FIG. 2: (a) Trajectories of vortices with the same circulation 
created at different radii but on the same polar angle, ob- 
tained by solving the GPE numerically; (b) Results obtained 
using the image method. 



equation for the position of the z-th vortex can be found: 




*,i = 1,2. (4) 

As observed in Figs. ^, ||, and the results from Eq. 
(m are in good agreement with the direct numerical sim- 



FIG. 3: Trajectories of a vortex-anti vortex pair. The vor- 
tices are initially created on the same radial distance but on 
opposing sides in the torus. We obtain the solid trajectories 
by solving the GPE numerically. The dashed trajectories are 
the computed solutions of the coupled differential equations 
(M). The period of oscillation for this setup is 225 ms. 
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FIG. 4: Trajectories of a vortex-anti vortex pair. The vortices 
are initially created on the same polar angle but on different 
radii. One vortex moves to the condensate boundary, the 
other one stays in the torus. We obtain the solid trajectories 
by solving the GPE numerically. The dashed trajectories are 
the computed solutions of the coupled differential equations 

(§)• 



ulations. Deviations from the numerical solutions of the 
GPE are basically due to deformations of the vortex core 
(see Sec. 0), since the image method does not take into 
account the internal structure of the vortex core. 
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Analytical results. Conformal mapping. 



by: 



In this section, we obtain the analytical expression of 
the vortex velocities explicitely, by employing a method 
which combines the previously presented image method 
and conformal mapping techniques |23[ . 



Analytical method 

We introduce a complex velocity potential p(z = x + 
iy) = (f>(x, y) + i"f(x, y), such that dp/dz — u — iw is the 
velocity v = (u,w). p(z) has to be analytic everywhere 
except for singularities, and therefore <p(x,y) and j(x,y) 
are both solutions to the Laplace equation. One of them 
can be chosen to be the real velocity potential : W<p = v. 
The curves <p = constant and 7 = constant intersect 
in right angles and correspond to the equipotential and 
stream lines. It can be easily seen, that 



p{z) = - — log (Z-Zj) 



(5) 



is the complex velocity potential of a single line vortex 
at position Zj, and consequently 



1 



v = 



2?r If- 



(6) 



In presence of rigid boundaries the velocity field is 
constrained to have a vanishing normal component at 
the boundaries and, therefore, since the flow is always 
in the direction of the stream lines, the boundaries cor- 
respond to curves of constant 7(0:, y). One can always 
arrange p(z) such that this constant is zero, so one has 
Im \p(z)] — for all z at the boundary, where Im indi- 
cates imaginary part. The potential p(z) is constructed 
from the Green's function of the corresponding boundary 
value problem, by superposing the contribution of the N 
vortices to the flow [M5l : 



N 

p( z ) = f{ z ) + ^Z^G{z lZj ) 

3 



(7) 



The term f(z) allows for some further flow due to mov- 
ing boundaries. The Green's function has the form 
G(z, Zj) = fund(z, Zj) + g(z, Zj), where fund(z, Zj) corre- 
sponds to the fundamental solution without any bound- 
ary (and also carries the divergence at the vortex posi- 
tions). According to (||): 



fund(z, Zj) 



2- 



log(z-Zj). 



(8) 



The function g{z,Zj) accounts for the constraints and 
it is smooth at the position of the vortices. It can be 
thought of as the complex potential induced by the im- 
ages. Therefore, the velocity of the j-th vortex is given 



Kj (uj - iwj) = Kj — (p{z) - fundO, 



dtp . dtp 
dyj dxj ' 



(9) 



where we have employed the, so called, Kirchhoff-Routh 
function E5I 



N 

i=i 



Im 



N 

2 

k=l 



G{zj,z k ) + -?-g{zj,Zj) 



(10) 



which is an integral of motion in the case of fixed bound- 
aries. Note, that contrary to the situation considered in 
electrostatics, there is no contribution from hxad(z, Zj) 
at the position of the vortex, since a vortex (possessing 
cylindrical symmetry) does not act on itself. 

For complicated geometries the velocity potential and 
the Kirchhoff-Routh function are difficult to calculate. In 
order to simplify one can map the complex plane which 
contains the complicated boundary value problem into 
another complex plane where the boundary geometry is 
simpler. A mapping F(z), which carries the z-plane into 
the C-plane, £ = F(z = x + iy) = s(x,y) + i t(x,y), is 
called conformal if F(z) is, apart from isolated points, 
an analytic function of z, the inverse mapping z = H(Q 
exists, and i?(£) is an analytic function of £• A confor- 
mal mapping maps continuous curves (in particular the 
boundaries) into continuous curves, and the solutions of 
the Laplace equation into solutions of the Laplace equa- 
tion, i.e. if the boundary value problem in the origi- 
nal plane is fulfilled, then it is also fulfilled in the new 
plane with the new boundaries beeing the mappings of 
the original ones. The transformation behaviour under 
the conformal mapping: 



K 2 



^ 47T 



dH 




~d( 


C=Cj 



(11) 



relates the Kirchhoff-Routh function of the old and new 
coordinates 



45 



This transformation law can be de- 
duced by the relation between the velocity of the j-th 
vortex in the z— plane and in the £— plane: 

dz-i 



dt 



dC 2 



C=Cj 



dt (dH/dOk^d 4tt 



((dff/dC)| 



(12) 



where Zj and Q denote the position of the j-th vortex in 
the corresponding coordinates and the bar denotes the 
complex conjugate. 
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One vortex in the torus 

Let us first calculate the velocity of one vortex in the 
torus. As already shown, the image configuration is 
slightly complicated. Let us now assume that the torus 
is centered at the origin of the slitted complex z-plane 
with radii R\ and i? 2 respectively. The mapping: 

Fi(z)=C = (2ir-6) + i(\n(r)-MR 1 )) 

= 2TT + ilogz-iliiR 1 (13) 

with 9 and r corresponding to the cylindrical coordinates 
of the z-plane, maps the annulus of the slitted z-plane 
onto a strip of width b — ln(i? 2 /i?i) and length 2ir in the 
( -plane, since we consider here only the principal branch 
of the logarithm. The inverse mapping is 



z = H 1 (Q = Ri expHC). 



(14) 



we get the Kirchhoff-Routh function of the original z- 
coordinates: 



lK*i) = — la 



sin 5- In — — 

llnf Ul 



+ ^la(n). 



(20) 



And from (|^) we obtain the vortex velocity: 



K I 7T 

cot 



4rln#^ lln#-""Vi?i 



(21) 



This analytical result is compared in Fig. [j] with the re- 
sults obtained from the direct summation provided by the 
image method. The agreement is, as observed, excellent. 



We approximate the velocity field in this strip of fi- 
nite length by the velocity field in a strip of infinite 
length. This approximation obviously becomes better as 
R2/R1 — * 1. Let us consider a second conformal mapping 

F 2 (C)=£ = expgc), ff 2 (£)=C=^log£ (15) 

which maps the infinite strip to the upper half plane and 
the boundaries of the infinite strip in the (-plane onto 
the real axis in the £-plane. We denote the position of 
the vortex in the (-plane as (1 and in the £-plane as 
£1. In the £-plane it is sufficient to consider one image 
vortex at position ^ to fulfill the boundary condition. 
The complex velocity potential in the £-plane is therefore 
given by 

P(0 = ~ log« - 6) + £ logtf - It). (16) 

Z7T Z7T 

The velocity of the vortex in the £-plane is: 



Ul - iW^ 



dt 



d 

in 



P(e) + gio g [£-£i] 



(17) 



According to ( jlzj ) we calculate the velocity of the vortex 
in the (-plane as: 



U i- iW i=Tb C ° t[ l tl] - 



(18) 



From this equation the Kirchhoff-Routh function in the 
C-plane can be obtained by solving (0): 



^( tl ) = J- In 



sin 



(19) 



Two vortices in the torus 

Let us now consider two vortices in the torus. We 
calculate first the velocity of the j-th vortex in an infinite 
strip of width b, containing two vortices. After applying 
again the mapping (|l5|), and following the same notation 
as in the previous section, the complex potential in the 
new coordinates is made up of the contributions of the 
two real vortices and two image vortices: 



P(0 



27T 8 U"6 



^logf^ 



(22) 



The velocity of the vortex 1 in the £-plane is given by: 

(23) 



tt( -rrrf iKl/2-IT iK2/2lT iK 2 /27T 

U{ - iWl — — 1 - - — 



& -s, a - £ 2 Si - 6 

According to ( |l2| ) we calculate the velocity of the vortex 
in the £-plane. After elemetary but tedious calculations 
we obtain: 



K\ fit 

77 cot 7 f 
4b \b 



11 
4b 



cosh (f (si - s 2 )) - COS (f (ti + t 2 )) 



sm 



:f (*a-*0) 



cosh (f (si - s 2 )) - cos (f (i a - *i)) 

sinh(f (gi-aa)) 
_i cosh (f ( Sl - s 2 )) - cos (f (t 2 - ti)) 

siah (I (s x - s 2 )) 



cosh (f (si - s 2 )) 



1 (*i+*a)). 

(24) 



where the ( denotes that this is the Kirchhoff-Routh func- 
tion for the (-plane. From the transformation law (111]) 



For the second vortex, we have just to interchange the 
indices 1 and 2. From (pi) we obtain the Kirchhoff-Routh 
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function: 



2 K 2 



4-7T 
47T 



4i ln ( sin l^" 



cos ( - (t\ - t 2 ) 



cosh 



COS 



(- b (t 1+h) ) 



cosh ( - (si - s 2 ) 



(25) 



and from (|llj) the Kirchhoff-Routh function for the orig- 
inal z-plane, corresponding to two vortices in the torus: 



ib = — In ( sin 

^ 47T 

J = l 



K\K2 
47T 

47T 



In 



In 



6 Ri 



fin H 

7r / rir 2 



2 K 2 

Air 

i=i 



Info) 



cosh ( — (0i 



o) 



ries as curves of constant Kirchhoff-Routh function (see 
Fig. §. 



ri,2 [IM] 

18r 

17 
16 
15 
14 
13 



cosh 



f- 



(26) 



Since the Kirchhoff-Routh function is for fixed bound- 
aries a constant of motion, we obtain from this equa- 
tion information about the vortex trajectories by com- 
puting the hyper surface on which the Kirchhoff func- 
tion is conserved. In particular we investigate the case of 
the vortex-antivortex pair, where the vortices are initially 
placed on the same radial coordinate. Due to symmetry 
considerations it is clear, that the vortices will always 
stay on the same radial coordinates r% — r 2 and with 
this knowledge it is easy to compute the vortex trajecto- 




0.1 0.2 0.3 0.4 0.5 0.6 



FIG. 5: Trajectories of a vortex-antivortex pair with vortices 
placed on the same radial coordinate. On the curves shown 
the Kirchhoff-Routh function is conserved. The vortices ap- 
proach, move in radial direction and separate again (see FIG. 
b|). The asymmetry of the image configuration is responsi- 
ble for the fact, that the trajectories are not symmetric with 
respect to reflection at r = 15/im. 

From the Kirchhoff-Routh function we get the velocity 
of the j-th vortex: 




(27) 
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with i,j — 1, 2, where i ^ j. In this formula one has to 
take into account, that the azimutal seperation between 
the vortices fulfills \6i — 6j\ < 7T. We have numerically 
solved the first-order differential equations resulting from 
this velocity field, and obtained the vortex trajectories, 
whose shapes are practically indistinguishable from the 
ones we obtain by numerically solving (0). 



III. HARMONIC TOROIDAL TRAP 

We discuss now the dynamics of vortices in a toroidal 
geometry with harmonic trapping in the radial direction: 



1 



Vtra P (r) = -zm (uP p (p - p ) 2 + u 2 z z 2 ) . 



(28) 



As before we consider a quasi-2D condensate in the 
Thomas-Fermi-Regime. In order to compare with the 
previous calculation for a box-like annulus, the trapping 
frequencies are adjusted such that the Thomas-Fermi 
radii correspond to an annular condensate with inner 
radius of 12 fim and outer radius of 18/im. The quali- 
tative dynamic of the vortices is unchanged from that in 
the box-like case, but the time scales are different. Due 
to the inhomogeneity of the trapping potential a vortex 
will no longer only move with the superfluid velocity. In 
addition it will have some velocity with respect to the 
superfluid, resulting from the small deviation from the 
cylindrical symmetric vortex state caused by the poten- 
tial. 

This velocity can be calculated by the method of 
matched asymptotic expansion |47| , Ell, consisting in 
evaluating the form of the solution of equation ([]]) in the 
vortex core and far away from it. The vortex velocity can 
be calculated by the requirement that the asymptotic be- 
havior of the two solutions has two coincide at a distance 
of the order of the healing length from the center of the 
vortex. The calculation yields |4q|: 




(29) 



where r*o denotes the position of the vortex and fyp is 
the Thomas-Fermi wave function. Since the deviation 
from the cylindrical symmetric vortex solution in the 
smooth potential is small, we neglect completely their 
effect on the velocity field induced by the vortex. There- 
fore, we can employ again the image method for deter- 
mining the velocity field fulfilling the boundary condi- 
tions. Therefore, the total velocity field results from 
adding to the field calculated in the previous section that 
of (^9|). This produces a very good agreement with our 
numerical results, as observed in Fig. |6|. 

The influence of the potential on the vortex velocity 
gives the possibility to control and manipulate the vortex 
dynamics by adiabatically changing the trapping poten- 
tial. For instance, it is possible to reverse the direction of 
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FIG. 6: Velocity of one vortex in the harmonic toroidal trap as 
a function of the radial coordinate. The solid curve shows the 
result of the numerical computation of the GPE, the dashed 
curve shows the combination of Eqs. ( pl| ) and (pfij). 



the vortex orbit in the annular trap (Fig. Q). Let us con- 
sider the above investigated boxlike toroidal trap, where 
we create one vortex. If one imposes a slowly increasing 
linear slope on the initially flat ground of the box the 
contribution of ( p9| ) eventually becomes more important 
then the boundary effects, and the vortex turns back- 
wards. 




FIG. 7: Trajectory of one vortex in a box like toroidal trap 
obtained by numerical solvation of the GPE. We add to the 
initially flat box like trap a linear, slowly increasing, slope. 
The vortex slows down as the slope increases and finally turns. 



The vortex-vortex interaction is also affected by the 
transversal harmonic potential, but the qualitative be- 
havior described in Sec. || is maintained J33| . 
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IV. DISSIPATION EFFECTS 

In this section, we discuss how dissipation affects the 
vortex dynamics in a toroidal geometry. Following a sim- 
ilar analysis as in Ref. |^8|, |4l[] we evaluate the free energy 
functional 



E 



£j|V*(r)f 



y(r)|^(r)| 2 



(30) 

for a condensate containing a vortex-anti-vortex pair 
with cylindrically symmetric vortex cores, for the vor- 
tex lines placed at the same radial coordinate r, but at 
different angles 9i,6 2 (see Fig. ||). Since the GPE de- 
scribes a Hamiltonian evolution and it does not include 
any dissipation, the trajectories of the vortices according 
to the GPE are curves of constant energy. However, if 
the vortices approach each other (or the boundaries) to 
a distance comparable to the healing lenghth, the vortex 
cores deform and the trajectories we obtain numerically 
slightly differ from these energy lines. The total energy 
increases with increasing distance between the vortices, 
or between a vortex and the torus boundaries, because 
the kinetic energy grows rapidly in this case. Therefore, 
we expect that in presence of dissipation, the vortices 
will tend to move towards the trap walls 0]. In addi- 
tion, the vortices will tend to reduce the distance between 
each other, and therefore, if the vortices do not collide 
with the walls, the apparent repulsion introduced by the 
toroidal geometry will be eventually overcome, and the 
vortices will eventually collide, and mutually annihilate, 
as in homogeneous condensates. 




geometry. The vortex trajectories in the presence of rigid 
boundaries have been obtained numerically by direct nu- 
merical computation of the corresponding 2D GPE. We 
have shown that these trajectories can be well described 
by closed analytical formulas, obtained from the com- 
bination of the image method and conformal mapping. 
Additionally, we have considered the case of harmonic 
transversal confinement, and showed that the numeri- 
cal results are in excellent agreement with the analyt- 
ical ones, once included the effects of the inhomoge- 
neous trapping potential calculated from the correspond- 
ing matched asymptotic expansion. 

We should note at this point, that the image method is 
applicable, when the vortices are spatially separated by a 
distance greater than their core size, which is of the order 
of the healing length. The latter statement also concerns 
the distance between the vortices and the boundaries. If 
this is the case, no deformation of the vortex core is pro- 
duced. However, if the vortices approach to distances 
comparable to the core size, the initially cylindrical sym- 
metric vortex cores experience evident deformations and 
therefore, the results provided by the image method be- 
gin to depart from the numerically obtained ones. Simi- 
larly, the dynamics of vortices in the quantum gas begin 
to differ significantly from the classical theory 1 44 , Et9j . 



The analysis of the vortex dynamics has been limited 
to the case of 2D traps. Lower dimensional condensates 
are currently actively investigated, and very recently the 
experimental observation of 2D (and even ID) condcn- 

the re- 



sates have been reported |37 



, |39J. However, 
suits of the present paper and Ref. |3j|, open interest- 
ing questions concerning the behavior of vortices in 3D 
toroidal condensates. Preliminar results concerning vor- 
tex lines parallel to the torus axis, confirm that in a 3D 
torus an essentially similar picture as that discussed for 
2D is obtained. However tiltings with respect to the torus 
axis, together with boundary effects, could lead to a sig- 
nificant bending of the vortex lines. This effect will af- 
fect the observability of the vortices, since the density 
minimum will depart from the direction of experimental 
observation. In addition, the bending will significantly 
distort the vortex dynamics, introducing an additional 
radial velocity to the vortex line. The vortex-vortex in- 
teraction is also expected to be significantly affected in 
the presence of a significant bending. The analysis of 
these effects, as well as those related with the propaga- 
tion and interaction of vortex rings in toroidal geometries 
will be the subject of further investigation. 



FIG. 8: Free energy of the condensate containing a vortex- 
antivortex-pair, depending on the positions of the vortices. 



V. CONCLUSIONS AND OUTLOOK 

In this paper, we have investigated the vortex dynam- 
ics and vortex-vortex interaction in an annular trapping 



Acknowledgments 

We acknowledge support from Deutsche Forschungs- 
gemeinschaft (SFB 407), TMR ERBXTCT-96-002, and 
ESF PESC BEC2000+. L. S. wishes to thank the Alexan- 
der von Humboldt Foundation, the Federal Ministry of 
Education and Research and the ZIP Programme of the 
German Government. Discussions with K.-A. Suominen 



10 



and J. P. Martikainen are acknowledged. 



M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. 
Wieman and E. A. Cornell, Science 269, 198 (1995); 
K.B. Davis, M. O. Mewes, M. R. Andrews, N. J. van 
Drutten, D. S. Durfee, D. M. Kurn and W. Ketterle, 
Phys. Rev. Lett. 75, 3969 (1995); 

C. C. Bradley, C. A. Sackett, and R. G. Hulet , Phys. 
Rev. Lett. 78, 985 (1997). 

D. G. Fried, T. C. Killian, L. Willmann, D. Landhuis, 
S. C. Moss, D. Kleppner, and T. J. Greytak, Phys. Rev. 
Lett. 81, 3811 (1998). 

S. L. Cornish, N. R. Claussen, J. L. Roberts, E. A. 
Cornell, and C. E. Wieman, Phys. Rev. Lett. 85, 1795 
(2000). 

A. Robert, O. Sirjean, A. Browaeys, J. Poupard, S. 
Nowak, D. Boiron, C. I. Westbrook, and A. Aspect, Sci- 
ence 292, 461 (2001); F. Pereira Dos Santos, J. Leonard, 
Junmin Wang, C. J. Barrelet, F. Perales, E. Rasel, C. 
S. Unnikrishnan, M. Leduc, and C. Cohen- Tannoudji, , 
Phys. Rev. Lett. 86, 3459 (2001). 

See e.g. E. M. Lifshitz and L. P. Pitaevskii, Statistical 
Physics, Part 2, (Pergamon Press, Oxford, U.K., 1980). 
O. M. Morago, S. A. Hopkins, J. Arlt, E. Hodby, G. 
Hechenblaikner, and C. J. Foot, Phys. Rev. Lett. 84, 
2056 (2000). 

R. Onofrio, C. Raman, J. M. Vogels, J. R. Abo-Shaeer, 

A. P. Chikkatur, and W. Ketterle, Phys. Rev. Lett. 85, 
2228 (2000) 

See e.g., R. J. Donnelly, Quantized Vortices In Helium II 
(Cambridge University Press, Cambridge, U.K., 1991). 
K.-P. Martzlin, W. Zhang, and E. M. Wright, Phys. Rev. 
Lett. 79, 4728 (1997); K.-P. Martzlin and W. Zhang, 
Phys. Rev. A 57, 3801 (1998); ibid, 4761 (1998). 
R. Dum, J. I. Cirac, M. Lewenstein, and P. Zoller, Phys. 
Rev. Lett. 80, 2972 (1998). 

B. Jackson, J. F. McCann, and C. S. Adams, Phys. Rev. 
Lett. 80, 3903 (1998). 

B. M. Caradoc-Davies, R. J. Ballagh, and K. Burnett, 
Phys. Rev. Lett. 83, 895 (1999). 

L. Dobrek, M. Gajda, M. Lewenstein, K. Sengstock, G. 
Birkl, and W. Ertmer, Phys. Rev. A 60, R3381 (1999). 
J. Williams and M. Holland, Nature 401, 568 (1999). 
D. L. Feder, C. W. Clark, and B. I. Schneider, Phys. Rev. 
A 61, 011601 (2000). 

K. W. Madison, F. Chevy, W. Wohlleben, and J. Dal- 
ibard, Phys. Rev. Lett. 84, 806 (2000); F. Chevy, K. W. 
Madison, J. Dalibard, Phys. Rev. Lett. 85, 2223 (2000) 
J. R. Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ket- 
terle, Science 292, 476 (2001). 

M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 
Hall, C. E. Wieman, and E. A. Cornell, Phys. Rev. Lett. 
83, 2498 (1999). 

B. P. Anderson, P. C. Haljan, C. A. Regal, D. L. Feder, 
L. A. Collins, C. W. Clark, and E. A. Cornell, Phys. Rev. 
Lett. 86, 2926 (2001). 

See e.g., A. L. Fetter and A. A. Svidzinsky, J. Phys.: 
Condens. Matter 13, R135 (2001), and references therein. 
G. K. Batchelor, An introduction to fluid mechanics, 
(Cambridge University Press, Cambridge, U.K., 1967) 



[24 
[25' 
[26 



P. O. Fedichev and G. V. Shlyapnikov, Phys. Rev. A 60, 
R1779 (1999). 

M. Guilleumas, and R. Graham, Phys. Rev. A 64, 033607 
(2001). 

E. M. Wright, J. Arlt, and K. Dholakia, Phys. Rev. A 
63, 013608 (20 00). 



[27] D. S. Rokhsar, |zond-mat/9709212 . 

[28] L. Salasnich, A. Parola, and L. Reatto, Phys. Rev. A 59, 
2990 (1999). 

[29] M. Benakli, S. Raghavan, A. Smerzi, S. Fantoni, and S. 

R. Shenoy, Europhys. Lett. 46, 275 (1999). 
[30] J. Tempere, J. T. Devreese, and E. R. I. Abraham, Phys. 

Rev. A 64, 023603 (2001). 
[31] J. Brand, and W. P. Reinhardt, J. Phys. B: At. Mol. Opt. 

Phys. 34, L113 (2001). 
[32] P. D. Drummond, A. Eleftheriou, K. Huang, and K. V. 

Kheruntsyan, Phys. Rev. A 63, 053602 (2001). 
[33] J. P. Martikainen, K. A. Suominen, L. Santos, T. Schulte, 

and A. Sanpera, Phys. Rev. A in press.; R. Dum, A. 

Sanpera, K.-A. Suominen, M. Brewczyk, M. Kus, K. 

Rzazewski and M. Lewenstein, Phys. Rev. Lett. 80, 3899 

(1998) 

[34] E. P. Gross, Nuovo Cimento 20, 454 (1961) 

[35] E. P. Gross, J. Math. Phys. 4, 195 (1963) 

[36] L. P. Pitaevskii, Zh. Eksp. Teor. Fiz. 40, 646 (1961) [Sov. 
Phys. JETP, 13, 451 (1961)] 

[37] A. I. Safonov, S. A. Vasilyev, I. S. Yasnikov, 1. 1. Lukashe- 
vich, and S. Jaakkola, Phys. Rev. Lett. 81, 4545 (1998). 

[38] A. Gorlitz, J.M. Vogels, A.E. Leanhardt, C. Raman, T.L. 
Gustavson, J.R. Abo-Shaeer, A. P. Chikkatur, S. Gupta, 
S. Ino uve, T.P. Rosenban d, D.E. Pritchard, and W. Ket- 
terle, |cond-mat/0104549| . 

[39] S. Burger, F. S. Cataliott i, C. Fort, P. Mad daloni, F. 
Minardi, and M. Inguscio, :ond-mat/0108037. 

[40] W. H. Press et al., Numerical Recipes, Cambridge Uni- 
versity Press, 1986 

[41] S. A. McGee and M. J. Holland, Phys. Rev. A 63, 043608 
(2001). 

[42] A. L. Fetter and A. A. Svidzinsky, J. Phys.: Condens. 

Matter 13, R135 (2001) 
[43] E. Lundh and P. Ao, Phys. Rev. A 61, 063612 (2000) 

and references therein. 
[44] A. L. Fetter, Phys. Rev. 151, 100 (1966) 
[45] P. G. Saffman, Vortex Dynamics, (Cambridge University 

Press, 1992). 

[46] Since we are dealing with a multiple connected region, the 
velocity field is not determined uniquely by the boundary 
condition and the specification of the vorticity field alone. 
This is reflected by the fact that the vortex centered at 
the origin can have arbitrary circulation without affecting 
the boundary condition. We choose the central vortex to 
have the circulation stated in the text, to prevent ring 
currents. Fixing in that way the circulation of the vortex 
at the origin, sets the so-called cyclic constant of the 
velocity field to zero. With this additional information 
the velocity field is determined uniquely jp3[ . 

[47] B. Y. Rubinstein and L. M. Pismen, Physica D 78, 1 
(1994). 



11 



[48] A. A. Svidzinsky and A. L. Fetter, Phys. Rev. A 62, [49] A. L. Fetter, Phys. Rev. 138, 429 (1965) 
063617 (2000). 



